Orbitally and Magnetically Induced Anisotropy in Iron-based Superconductors 
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Recent experimental developments in the iron pnictides have unambiguously demonstrated the 
existence of in-plane electronic anisotropy in the absence of the long-range magnetic order. Such 
anisotropy can arise from orbital ordering, which is described by an energy splitting between the 
two otherwise degenerate d^z and dyz orbitals. By including this phenomenological orbital order 
into a five-orbital Hubbard model, we obtain the mean-field solutions where the magnetic order is 
determined self-consistently. Despite sensitivity of the resulting states to the input parameters, we 
find that a weak orbital order that places the dy^ orbital slightly higher in energy than the d^z orbital, 
combined with intermediate on-site interactions, produces band dispersions that are compatible with 
the photoemission results. In this regime, the stripe antiferromagnetic order is further stabilized 
and the resistivity displays the observed anisotropy. We also calculate the optical conductivity and 
show that it agrees with the temperature evolution of the anisotropy seen experimentally. 

PACS numbers: 74.70.Xa, 75.25.Dk, 71.10.Fd 
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I. INTRODUCTION 

In-plane anisotropy plays a prominent role in iron- 
based superconductors. For example, inelastic neutron 
scattering (INS) first established that the exchange con- 
stant along the ferromagnetic y direction, Jiy, is not only 
much smaller than the one in the antiferromagnetic x 
direction, Ji^,, but is actually slightly negative^i sug- 
gesting that an unfrustrated spin model underlies the 
magnetism. Angle-resolved photoemission spectroscopy 
(ARPES) also observed a dramatic orbital-dependent 
Fermi-surface reconstruction upon the magnetostruc- 
tural phase transitionj^i^ However, due to the fact that 
the crystals used in such experiments spontaneously form 
dense domains, the signals from the two diagonal phases 
were mixed in these early experiments. Hence, it was cru- 
cial that scanning tunneling microscopy (STM) detected 
a quasi-one-dimensional interference pattern,^ thus con- 
firming that the anisotropy arises entirely from a sin- 
gle domain. Helping to complete the story were the 
breakthrough experiments in which a uniaxial stress was 
applied to almost fully detwin the crystals <^ This tech- 
nique revealed an in-plane anisotropy in both the do^^— 
and ao^"— conductivities. Consistent with this transport 
anisotropy are the band dispersions derived from ARPES 
measurements on detwinned samples.— 

The observed anisotropies (for a review, see Ref. fisl ) 
have distinct origins in terms of the local-itinerant elec- 
tron dichotomyfiir— On the one hand, INS measures 
the spin excitations that arise predominantly from lo- 
cal moments. Although a structural transition pre- 
cedes or coincides with the onset of the antiferromag- 
netic order, its magnitude is too small to account 
for the strong anisotropy. In this context, several 
theoriesii"— adopted a Kugel-Khomskii^i model with 
orbital-dependent superexchanges and proposed that 
ferro-orbital order stabilizes the (tt, 0) magnetic order, 
leading naturally to anisotropic spin couplings. In an 
alternative approach j^^i^^ the local moments, which are 



governed by an isotropic J1-J2 Heisenberg model, are 
coupled with the itinerant electrons of the two degen- 
erate dxz and dyz orbitals by the Hund exchange. From 
the double-exchange mechanism, a ferro-orbital order in 
the itinerant band leads to different corrections to the 
spin-wave excitations along two diagonal directions, re- 
sulting in an anisotropic spin-only model. 

On the other hand, the detected anisotropy could sim- 
ply signify that the itinerant electrons also undergo a 
phase transition that breaks the C4 rotational symme- 
try, which is no surprise since the symmetry is already 
broken by the underlying (tt, 0) antiferromagnetism. Al- 
though this electronic anisotropy does persist above the 
Neel temperature and the structural transition tem- 
perature Ts in some experiments, ^"^^ it should be noted 
that the applied uniaxial stress unfortunately turns the 
sharp transition into a broad crossover. In this sense, re- 
cent INS measurements^ that observed a spin nematic 
order in the paramagnetic phase have taught us some- 
thing new - the magnetic anisotropy, and the possible 
electronic anisotropy, can exist in the absence of the (tt, 0) 
antiferromagnetism and even the structural distortion. 
This "nematic" order is indeed confirmed by the band 
splitting above T5 observed by the latest ARPES exper- 
iments on twinned samples,^^^, and recently hinted at by 
the zero-bias conductance enhancement in point-contact 
spectroscopy.^^ More convincing evidence comes from the 
anisotropic vortex core structure in the superconducting 
FeSe,— which exhibits no static spin order. 

One theoretical proposal that is in accord with some 
of the current experimental observations is the J1-J2 
Heisenberg model, in which an Ising order can occur 
without long-range magnetic order ;^^!^^ according to the 
"order from disorder without order" mechanism;^^' How- 
ever, the unfrustrated magnetism observed by INS^i^^ 
suggests that the correct physical model lies elsewhere. 
Another possible candidate is a Kugel-Khomskii spin- 
orbital model^— that might support a phase which pos- 
sesses some orbital order that breaks the C4 rotational 
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symmetry, but no long-range magnetic order. Be- 
sides these strong-coupling approaches, there are other 
theoretical proposals that attribute the nematic or- 
der to interaction-induced instability of the itinerant 
electrons. ^^"^^ 

Our particular view, which we have enunciated 
elsewhere,—!^ is that orbital order (rather than any in- 
trinsic tendency of the electrons to orient themselves 
inhomogeneously) lies at the heart of the anisotropics 
and the observed structural and magnetic transitions. 
Certainly, it is now commoD^i. to refer to any phase 
that breaks C4 symmetry as an electron nematic, even 
when such a phase is lattice or orbitally induced. 
This view has been scrutinized sharply through re- 
cent experiments^ on the the cuprate superconductor 
Bi2Sr2Can_iCu„02nH-4+a;, which indicate that many of 
the anisotropics that have been viewed as evidence for 
an electron nematic phase actually originate from lattice 
effects in the BiO layer. This experiment certainly indi- 
cates that caution rather than a rush to nematize is war- 
ranted. For the pnictides, the origin of the anisotropics 
will remain open in the absence of a clear experiment that 
is able to discern their efhcient cause. Consequently, our 
usage of the term nematic here strictly entails the orbital 
order which necessarily breaks the equivalence between 
the canonical x and y axes. 

In this paper, we will not focus on the origin of this 
"nematic" order (in the sense defined above), but rather 
investigate its experimental consequences by modeling it 
phenomenologically as an energy-splitting term between 
the Fe dxz and dyz orbitals. This type of orbital ne- 
matic order has been studied previously,^- but only in 
the context of an orbital- independent magnetic order, 
which is insufficient to capture the complicated electronic 
structure. To this end, we will start with the multi- 
orbital model that explicitly includes this orbital ne- 
matic order and solve the mean-field Hamiltonian using 
a self-consistent procedure. In this approach, the mag- 
netic order on different orbitals will be determined more 
realistically by the band structure and the interaction 
strength. Moreover, we will be able to address how the 
orbital and magnetic orders interplay with one another. 
To reach agreement with photoemission experiments 
we find that the d^z orbital is placed lower in energy 
than dyz, and that this orbital order strengthens the 
stripe antiferromagnetism. The orbital and magnetic or- 
der together reconstruct the band structure and result 
in the anisotropy at both the low-temperature antiferro- 
magnetic and the high-temperature paramagnetic phase. 
These findings suggest that orbital order plays the central 
role in the electronic structure of the iron-based super- 
conductors. 

The paper is organized as follows. The general for- 
malism of the multi-orbital model is described in Sec. |lll 
We introduce the orbital nematic order in Sec. IIIII and 
calculate its experimental consequences that are relevant 
for anisotropy in the paramagnetic phase. Section IIVI 
presents the complicated electronic structure arising from 



both the orbital and magnetic orders. It is shown that 
both of them are indispensable components that lead to 
anisotropics observed by many experiments in the mag- 
netically ordered state. The implications of our findings 
are discussed in the last section. 



II. MULTI-ORBITAL MODEL 

The multi-orbital Hamiltonian we start with is usu- 
ally defined within an extended Brillouin zone that only 
contains one Fe atom per unit cell. The kinetic-energy 
contribution is written as 

where c\ creates an electron of momentum k on or- 

bital a with spin /i (/i =tj i)- We note here that fc is not 
the crystal momentum k defined by the translation op- 
erator % of the Fe lattice unit vector {i — x, y, z), i.e., 
Ti\k) = e^^^ \k). Rather, it is the operator VzTi, where Vz 
is the reflection operator in the z direction, instead of 7i, 
that leaves the Fe-As lattice invariant. Thus, k actually 
labels the eigenstates oiVzTi, i.e., VzTi\k) — e'^^'lfc), and 
yields an unambiguous way to unfold the real Brillouin 
zone with two Fe atoms per unit cell. This important dis- 
tinction between fe and k was previously discussed in de- 
tail in Ref.lsrl In principle, we need to fold back the band 
dispersions that are obtained by the diagonalization of 
T~Lk (HJ, and fit them to the local-density-approximation 
(LDA) calculations by tuning the tight-binding hopping 
parameters. This has been done in many studies, using 
two^-iS three,.^i^ four^ or five^^i^ Fe d orbitals. In 
order to make our calculations more realistic, we only 
focus on the five-orbital model, particularly the one of 
Ref. 113, which is based on a three-dimensional fitting to 
the LDA band structures of BaFe2As2, the material on 
which most of the experiments are performed. 

It is helpful for us to return to real space where T-Lk 
(P) takes the form Hk = J2itf 4atj.^jfit^^ where t"f is 
the hopping amplitude, with ?, j denoting the index of 
the site. As discussed, the operator PzTi leaves Hk in- 
variant. Under Vz%, we have Cia^ XaCi+ei,ati, where 
= 1 for a = dxy, dx2_y2, or d3^2_^2, and Xa = ~1 
for a — dxz or dyz- Thus, it is required that t"^^ = t"fj 
for XaXfs = 1, whereas t"/ = e'^-'^'-t'^!^- for XaXp = -1, 
where K = (7r,7r,7r) and i"^^- only depends on — rj. 
We immediately noticed that in the crystal momentum 
space, the electron operators of the dxz and dyz orbitals 
at k are coupled with those of dxy, dx2^y2, and d-^z^_j.2 at 
k-\- K . As a result, we define the pseudocrystal momen- 
tum k as follows: c^,^^ = Cka^i for a = dxz or dyz, and 
Cfca^ = Ck+K,afi for a ^ dxy, dx2^y2, or d-i^2_^2M^^^ 
It is in this pseudocrystal momentum k space that the 
kinetic energy T-Lk takes the diagonal form of %k 
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The distinction between k and k has immediate con- 
sequences for the interpretation of the ARPES measure- 
ments. The momentum probed by ARPES is not fc, but 
fc. Because the crystal momentum k is not a good quan- 
tum number of Hk © , ARPES detects both bands with 
momentum k and fc -I- K, corresponding to folding the 
Brillouin zone with one Fe atom per unit cell by a wave 
vector K — (tt, tt, tt). This observation is consistent with 
the fact that the real unit cell including the As atoms 
consists of two Fe atoms. However, the relative inten- 
sity of each band measured by ARPES depends on the 
strengths of hybridizations between fc and k + K. 

Formally, the Green function in terms of fe is defined 

as 

GMk,r)^-{TrCj^^^{T)cl_^^{0)), (2) 

with its Fourier transform satisfying Q^^{k,uj) — ujI — 
^(fc). We can simply write down the spectral function 



A{k,uj) = --'^ irg{k,uj + i5) 

TT 



(3) 



However, what is really measured by ARPES is not 
A{k,uj), but 



A{k,uj) 



1 



E 



Qaa{k + K,UJ + id) 



(4) 



It is this function on which we will focus. 

In Figs. [Tfa) and (b), we plot A{k,uj — 0), depicting 
the Fermi surfaces in the plane oi kz — and — n, re- 
spectively. Those of A{k,ui = 0) are shown in Figs.[TJc) 
and (d). The frequency w is defined with respect to the 
chemical potential n, which is determined by the elec- 
tron filling level n — 6, dictated by the Fe valence in 
the parent compounds. In addition, all the energy scales 
have an implicit unit of eV throughout this paper. By 
comparing these graphs, we find that A{k,uj) displays a 
more complicated structure due to the folding of the Bril- 
louin zone. From now on, we will only focus on A{k,Lo), 
which is measured by ARPES in our direct comparison 
with experiments. Furthermore, in Fig. [IJe), A{k,uj) is 
plotted along the line X-T-X, with X = (7r,0,0) and 
r = (0, 0, 0), to represent the band dispersions probed by 
unpolarized ARPES. Experimentally, the orbital charac- 
ter can be investigated by tuning the polarization of the 
incoming light. For example, using a polarization per- 
pendicular to the incident plane only selects those or- 
bitals that are odd under Vy {dyz and d^y) along the 
line of X-T-X. We plot these orbital-polarized disper- 
sions in Fig.[IJf), which shows qualitative agreement with 
experiments.— Since the C4 rotational symmetry is re- 
spected, the dispersions are exactly the same along the 
crystal momentum line Y-T-Y, where Y = (0, tt, 0), with 
no splitting between the bands at X and Y . 
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FIG. 1. (Color online) Fermi surfaces in terms of (a), (b) the 
pseudocrystal momentum fc and (c), (d) the crystal momen- 
tum fc. We plot the Fermi surfaces in the xy plane and set 
the z component in each figure to (a) fcz = 0, (b) fcz = tt, (c) 
kz = 0, and (d) kz = tt. (e) Unpolarized dispersions and (f) 
polarized dispersions along the crystal momentum line of X- 
V-X. In (f), only the components of the dyz and dxy orbitals 
are shown. We represent the values of the spectral function 
A{k,Lo) by the color scale, which is used consistently for all of 
the figures in this paper. 



III. ORBITAL NEMATIC ORDER 



As discussed in Sec. U electronic anisotropy has been 
confirmed by recent experiments^^iSi to persist above the 
magnetostructural transition. Our take on this is that 
this effect is due entirely to orbital ordering. To test out 
this hypothesis, we introduce the orbital nematic order 
as an energy-splitting term between the dxz and dyz or- 
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FIG. 2. (Color online) Fermi surfaces in the plane for (a) 
kz = and (b) kz = n. Polarized dispersions along the line of 
(c) X-F-X and (d) Y-T-Y. We have set the orbital nematic 
order parameter A = —0.08, explicitly breaking C4 symmetry. 

bitals, 

where = ±A for a = d^z and dy^, respectively, and 
Aq = for the other three orbitals. In principle, all five 
orbitals should be involved in this nematic order. But we 
will only consider the d^z and dyz orbitals due to their 
spatial anisotropy, whereas the other three orbitals are 
dropped because they are C4 rotationally symmetric. It 
needs to be emphasized that "Hat © represents an elec- 
tron nematic order which occurs without the onset of the 
long-range stripe antiferromagnetism. We find that in or- 
der to produce results that are consistent with ARPES 
measurements the orbital nematic order parameter A 
is required to have a small negative value, which leads to 
a higher energy of dyz relative to d^z- For the purpose of 
illustration, we choose A = —0.08, and plot the resulting 
Fermi surfaces and polarized dispersions in Fig. [51 

From Figs.llja) and (b), we notice that the hole pock- 
ets at the zone center are deformed into elliptical shapes 
that are elongated along the x and y direction in the 
plane of fc^ = and fc^ = tt, respectively. The apparent 
breaking of C4 symmetry is indeed a direct consequence 
of 'Hjv ® that is explicitly introduced by hand. How- 



ever, this orbital order does have a physically reasonable 
origin, which can be explained as follows. From a weak- 
coupling point of view, the observed magnetic order with 
wave vector q = (tt, 0, tt) arises from the nesting insta- 
bility between the hole pockets centered at A; = (0, 0, kz) 
and the electron pockets at A; = {TT,0,kz + tt). But the 
nesting is not perfect because the hole pockets are more 
circular, whereas the electron pockets are more ellipti- 
cal [see Figs. [T{a) and (b)]. By the inclusion oi Hn ^ 
with a small negative A, the hole pockets are deformed 
into ellipses, whereas the electron pockets are relatively 
less affected, thus resulting in a better nesting condition 
between the two by the wave vector q [see Figs. [2][a) 
and (b)]. Consequently, this type of orbital order will 
naturally arise in the system and provide further stabi- 
lization of the antiferromagnetism. Our result is con- 
sistent with the Pomeranchuk instability from the func- 
tional renormalization-group studiesi^ 

The polarized dispersions along the x and y direction 
are displayed in Figs.ElJc) and (d), respectively. It needs 
to be emphasized that only the spectral functions of dyz 
and dxy orbitals are shown in Fig.[2fc), whereas we only 
plot those of dxz and dxy in Fig. [2][d) , which are exactly 
what are measured by the polarized ARPES setup.l- In- 
deed, a small negative A, which lifts dyz higher than d^z 
in energy, produces splitting between the bands at X and 
y, in agreement with experimental observations. 

To further examine the validity of the orbital order 
term Hn ([5]) , we turn to its experimental consequences in 
transport measurements by calculating the optical con- 
ductivity (Ta;a;(a;) and ayy{uj) along the x and y direc- 
tion, respectively. When A = 0, axx and ayy are com- 
pletely equivalent [Fig. (JJa)], preserving the C4 symme- 
try. We see a sharp Drude peak accompanied by some 
high-energy features due to interband transitions. If a 
small negative A = —0.08 is turned on, a{uj) along 
two diagonal directions become distinct and a low-energy 
subpeak emerges around loi ~ 0.2, where ayy has a value 
larger than axx , as depicted in Fig. ^h) . We also notice 
that around the high-energy peak at i02 ~ 0.7, axx is 
dominant instead. These results are in good agreement 
with experiments^ performed above Tpf in the detwinned 
samples. 

We also compare the Drude weight of axx and ayy 
and plot their ratio as a function of A. As shown in 
the inset of Fig. ^a) , when the energy splitting A gets 
larger, axx{0) / cryy{0) almost monotonically decreases, re- 
sulting in a larger Drude weight along the y direction. 
If we naively assume proportionality between the Drude 
weight and the dc conductivity, this result contradicts 
experimental findings in which the antiferromagnetic x 
direction always has a higher conductivity^ than the fer- 
romagnetic y direction. However, as pointed out by op- 
tical measurements)^ the Drude weight does obtain a 
higher value along the y direction, and it is the scatter- 
ing rate difference that dictates the higher conductivity 
along the x direction. In this regard, to obtain the cor- 
rect anisotropy of the dc conductivity, we need to take 
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FIG. 3. (Color online) Optical conductivity a^x and ayy as 
a function of the frequency uj for the orbital nematic order 
parameter (a) A = and (b) A = —0.08. The inset of (a) 
displays the ratio of the Drude weight, o-xx{0)/a-yy{0), as a 
function of A. (t(ij;) is plotted in an arbitrary unit, which is 
kept the same in this paper. The two arrows in (b) denote the 
two characteristic frequencies cji « 0.2 and UJ2 ~ 0.7, where 
a{uj) exhibits a peak structure. 



IV. MAGNETIC ORDER 

Now we set out to include magnetic order. At the 
mean-field level, magnetic order can be described by 



■Hai 



al3 i 



.iqri z T „ 



(6) 



for a multi-orbital system, where we have chosen the or- 
dering wave vector q = (tt, 0, tt) and the spins are as- 
sumed to point along the direction, with cr^ being 
the Pauli matrix. For a five-orbital system, the order 
parameters Map form a 5 x 5 Hermitian matrix, with 
a total of 25 independent real variables. However, as 
discussed in Ref. [S^, the (tt, 0, tt) magnetically ordered 
state is invariant under the inversion X around any Fe 
site, reflections Vx, Vy, and Pz along the x, y, and z 
directions, respectively, and effective time reversal TTZ' , 
which is a combination of time reversal and spin rever- 
sal, TTZ' = TTZ o STZ. Under these symmetries in our 
five-orbital system, only six parameters acquire nonzero 
real values: Maa {oc = dxz,dyz,dxy,dx^-y2,d^z^_r2) and 
Mai3 = Mfja (a = dx2-y2, j5 = d^z^-r^)- Heuce, mag- 
netic order obtains almost exclusively within the same 
orbitals, with the only exception being the orbital-off- 
diagonal term between the d^^-y^ and d^z'^-r^ orbitals. 

By following these discussions, we consider the on-site 
interaction of a multi-orbital Hubbard model. 



J 



J I -j- 



J' 



Et t 



(7) 



where hi 



into account some anisotropic scattering mechanism, for 
example, a short-range (tt, 0,7r) magnetic order which is 
supported by a recent INS experiment.^'* Indeed, this cor- 
rect anisotropy of the scattering rate is obtained by a 
recent theory^ that considers scattering by anisotropic 
spin fiuctuations in the paramagnetic phase. It should 
also be mentioned that our result of the Drude weight 
agrees with earlier studies^ based on the same form of 
the orbital order, but using a different five-orbital tight- 
binding model4i 



In conclusion, our study shows that orbital order plays 
an important role in modifying the electronic structure 
as well as the ratio of the Drude weights, (Txx{^)/ <yyy{^)i 
above the onset temperature for long-range magnetic or- 
der. In this way, we have separated the physics that 
results from magnetic as opposed to orbital order. 



and ,] = ,]' 



We stiff assume U = V + 2J 
which is not necessarily valid as the orbitals 
used here only share the same symmetry but do not have 
exactly the identical form of the atomic d orbitals. We 
make use of the standard mean-field decoupling. 



(8) 



where /i = ±1 for up and down spins, respectively. As 
shown by LDA calculations;^ the Fermi surfaces are 
mostly composed of the t2g orbitals [dxz, dyz, and dxy). 
Thus we can safely ignore the orbital-off-diagonal mag- 
netic order between dx2-y2 and d^z^-r^i Sind use the 
above orbital-diagonal decoupling which captures five of 
the total of six nonzero mean-field antiferromagnetic or- 
der parameters. 

By contrast, we point out that in Ref. [4^, the orbitals 
are defined along the axes of the original unit cell, X and 
Y, which are rotated by 45° from the x and y axes of the 
Fe lattice. Applying mean-field decoupling of Eq. ([5]) on 
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this model will only take account of four order parame- 
ters. Symmetry considerations impose that the orbital- 
diagonal elements Maa of the dxz and dyz orbitals are 
equal to each other. However, the off-diagonal element 
Mai3 = Mj3a (a — dxz, /3 = dYz) can acquire nonzero 
values, but will not be captured by the mean-field the- 
ory. Hence, we will use the model'*^ where the orbitals 
are defined along the Fe-Fe bond, and expect better re- 
sults compared to earlier studies^i^ based on Ref. [H 
under orbital-diagonal mean-field decoupling. 

Straightforward calculation yields the mean-field inter- 
action term 



where 




and the constant 

+ ^'Yrnamj3. (12) 

Note that in Hj (0) , k can be simply replaced by k with- 
out changing the form of the equation. The full Hamil- 
tonian Ji = Hk + Hn + is quadratic in electron 
operators c^^^ , and can be solved with order parameters 
ria and rria being determined self-consistently. 

So far there is still no systematic mean-field study 
of this three-dimensional tight-binding model^^ specially 
constructed for BaFe2As2. Thus, as the first step, we 
need to search for appropriate values of U and J that are 
consistent with experimental observations. At the out- 
set, we set the orbital nematic order A to zero to sim- 
plify our calculations. Generally, a nonzero A produces 
further modifications, but the physically relevant regime 
of U and J is not greatly affected by the choice of A. 
In Fig. m we plot the total staggered magnetic moment 
TO = TOq as a function of Coulomb repulsion U for 
various Hund's couplings J . It is found that there exists 
a metallic phase with antiferromagnetic order at inter- 
mediate Coulomb repulsion U . Furthermore the ratio of 
J and U also needs to take intermediate values. This 
requirement of U and J is qualitatively consistent with 
earlier studie a^^'^^ based on other five-orbital models. 

For a typical value of U = 1.08 and J = 0.20C/, we 
plot the polarized dispersions along the two diagonal di- 
rections in Figs. [SJa) and (b), respectively. Clearly, the 
magnetic order opens up a gap close to Fermi energy and 
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FIG. 4. (Color online) The total staggered magnetic moment 
m as a function of the Coulomb repulsion U for different 
Hund's exchanges J. We set the energy splitting A = in 
order to find the regime of parameters that are of interest in 
the context of iron-based superconductors. 




Y r Y 

FIG. 5. (Color online) Polarized dispersions along (a) X-V-X 
and (b) Y-V-Y, for U = 1.08, J = 0.2017, and A = 0. The 
corresponding magnetic moment m = 0.46. 



significantly modifies the original band structure shown 
in Fig. [1] However, this (tt, 0, vr) antiferromagnetism, al- 
though breaking the C4 symmetry, produces only a small 
splitting between the bands at X and Y . In fact, the 
band energy at Y is even slightly higher than that at 
X. From our calculation, using other values of U and J, 
this near degeneracy of the bands at X and Y is quite 
robust, which contradicts ARPES measurements As 
will be seen, this problem can be alleviated by includ- 
ing orbital ordering. Hence, orbital order Hn © is still 
present in the antiferromagnetic phase and has important 
experimental consequences. 

In order to simplify our discussions, we set J = 0.20[/ 
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FIG. 6. (Color online) The total staggered magnetic moment 
as a function of the Coulomb repulsion U and orbital nematic 
order parameter A. We set Hund's exchange J = 0.20f/. 



and focus in the regime [/ e [1.0,1.2] and A G [—0.1,0.1], 
where significant agreement with experiment can be 
found. We first investigate variations of the total stag- 
gered magnetic moment m for different U and A. As 
discussed in Sec. IIIIl a negative A induces better nest- 
ing between hole and electron pockets and stabilizes the 
(tt, 0, tt) magnetic order, thus leading to an increase of 
m, as shown in Fig. [SI In contrast, the magnetic mo- 
ment m decreases when A takes larger positive values. 
In fact, a positive A places the dxz orbital higher in en- 
ergy compared to dyz , and favors the antiferromagnetism 
with q — (0,7r,7r) instead. This interesting interplay 
between the orbital and magnetic order was also noted 
recently^'' based on LDA results and a Ginzburg-Landau 
phenomenological theory. 

To further illustrate the role of orbital order in the 
reconstruction of the electronic structure in the magneti- 
cally ordered phase, we plot the Fermi surfaces and polar- 
ized dispersions in Fig. [T] We use the set of parameters, 
U = 1.08, J = 0.20t7, and A = -0.08, to be consis- 
tent with previous figures. From Figs. [TJc) and (d), the 
desired band splitting between X and Y is successfully 
produced due to a negative A. Furthermore, there are 
multiple band crossings at the Fermi energy along X-T- 
X, whereas only a single crossing occurs along Y-T-Y, 
also in agreement with ARPES findingsji2, We further 
point out that the two small Fermi surfaces adjacent to 
the large hole pocket in Fig.[7|b) are actually Dirac cones, 
which have been predicted theoretically^^ and confirmed 
experimentally] ^^'^"^ Note, however, that the existence of 
Dirac cones relies on the degeneracy between the d^z and 
dyz orbitals/^^ But the small orbital order used here is 
not enough to annihilate such Dirac features. Finally, 
the two Fermi-surface segments close to the zone center 
in Fig. EJa) are mostly aligned along the ferromagnetic 
y direction and in principle can produce the quasi-one- 
dimensional interference pattern observed in STMi^ 
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FIG. 7. (Color online) Fermi surfaces in the plane of (a) 
kz — and (b) kz = n. Polarized dispersions along the line 
of (c) X-F-X and (d) Y-T-Y. The parameters used here are 
U = 1.08, J = 0.20f/, and A = -0.08, with the total mean- 
field staggered moment m — 0.62. 
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FIG. 8. (Color online) Polarized dispersions along (a) X-T-X 
and (b) Y-r-Y, for U = 1.13, J = 0.2017, and A = 0. The 
corresponding magnetic moment m — 0.61, which is close to 
m = 0.62 of Fig. [71 
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FIG. 9. (Color online) (a) Optical conductivity a^xiijj) and 
ayy{uj) calculated using U = 1.08, J = 0.20(7, and A = -0.08. 
They are plotted in the same unit as Fig. [31 (b) Anisotropy 
of Drude weight, axx{0) /ayy{0), as a function of U and A. 
The black line separates the regions where a^xiO) > o-yy{0) 
and where axx{0) < ayy{0). We choose J = 0.2017 here. 



tivity shows two peaks at frequency uji and UJ2, which 
have similar values as we defined previously. Compared 
to the results without the magnetic order (see Fig. [3]), the 
magnitude of the peak at uji increased significantly. So 
the onset of this low-energy peak at wi is mostly due to 
antiferromagnetism, which opens up a gap at the chem- 
ical potential and shifts part of the Drude weight to the 
region around lui. On the other hand, the high-energy 
peak at uj2 comes from the interband transitions that 
are already present in the original multi-orbital model. 
However, we cannot distinguish along which direction 
the optical conductivity dominates since there are multi- 
ple crossings (as a function of frequency) between the two 
curves of cr^x and ayy. But, naively, a^x does have higher 
values in most of the regions around wi , which matches 
the experimental results^ at the low-temperature antifer- 
romagnetic phase. Also of note is the region around uji 
where a^x displays a single peak, whereas ayy acquires a 
double-peak structure, which is also reproduced in a re- 
cent study^ using the combination of density functional 
theory and dynamical mean- field theory (DFT-I-DMFT). 
Anisotropy of the Drude weight is presented as a func- 
tion of U and A in Fig. Hljb). If we assume an isotropic 
scattering rate, the dc conductivity anisotropy has the 
correct directionality in a very large region of parameter 
space. However, considering the (tt, 0, tt) antiferromag- 
netism, we should expect a larger scattering rate along 
the X direction, where the spins are antiparallel to their 
neighbors. This should further restrict the parameter 
space of interest. Nonetheless, we find that in a rea- 
sonable range of parameters, both orbital and magnetic 
order underlie the resistivity anisotropy in a multi-orbital 
model of the pnictides. 



However, attentive readers may notice that for Fig. [71 
the corresponding magnetic moment m — 0.62, which is 
larger than m = 0.46 of Fig. [51 Hence, it is possible that 
the splitting between the bands at X and Y is caused 
by the stronger magnetic order in Fig. [71 In order to 
confirm the splitting is entirely due to the orbital order 
term Hn ([S]), we plot the polarized dispersions in Fig. [H 
for U = 1.13, J = 0.20C/, and A = 0. These parameters 
produce a staggered magnetic moment m = 0.61, which 
is close to the value of m = 0.62 in Fig. [71 However, 
no spUtting between the bands at X and Y is generated 
in Fig. [SI where the band energy at X is even slightly 
lower than that at Y. As we have pointed out earlier, 
this near degeneracy between the band energy at X and 
Y survives for other values of U and J as well, as long as 
the orbital nematic order parameter A = 0. Therefore 
we have shown that it is the orbital order Hn (0), rather 
than the stripe antiferromagnetism, that is responsible 
for the band splitting at X and Y observed in ARPESJ^ 

We end this section by discussing the optical spec- 
tra in the presence of both orbital and magnetic order. 
The same set of parameters as those of Fig. [71 are used 
here. From Fig. ^a), we see that the optical conduc- 



V. FINAL REMARKS 

In summary, we have solved a three-dimensional five- 
orbital model using the mean-field approximation in the 
presence of both orbital and magnetic order. We mod- 
eled the orbital order phenomenologically as an energy 
splitting term between the two otherwise degenerate d^z 
and dyz orbitals, and find that it is solely responsible for 
the splitting between the bands at X and Y observed 
by polarized ARPES.^^ This orbital nematic order also 
causes the Drude weight difference between the two diag- 
onal directions for the paramagnetic phase, in agreement 
with the analysis based on optical conductivity;^ On the 
other hand, the magnetic order sets in at a lower temper- 
ature. It opens up a gap at the Fermi energy and shifts 
part of the Drude weight to high energies, leading to the 
emergence of a peak structure at uji ~ 0.2 in the opti- 
cal conductivity. However, magnetic order alone cannot 
explain the anisotropic dispersions, especially the band 
splitting between X and Y. Thus both orbital and mag- 
netic orders are present in the low-temperature antiferro- 
magnetic phase, and together induce various anisotropics 
seen by different experimental techniques. This result 
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is in contrast with earlier studieS j^^'^^i^ which claimed 
that the magnetic order is the driving mechanism for the 
anisotropy. Thus, our study establishes the importance 
of orbital order, irrespective of whether long-range mag- 
netic order is present or not. 

Besides orbital order, the anisotropy observed above 
the magnetostructural transition can come from a differ- 
ent mechanism, for example, static short-range magnetic 
order. Its existence is indeed supported by recent INS 
experiments.-^^ In principle, we can model this by a prob- 
ability distribution of the ordering wave vector peaked 
a.t q = (tt, 0, TT)r^ and investigate its modification to the 
electronic structure. It is expected that the magnitude 
of the subpeak at uji in Fig.[3ljb) will increase due to this 
short-range magnetic order, thus yielding better agree- 
ment with experiments.^ Nevertheless, the orbital order 
should always be present until the band splitting van- 
ishes. 

Of course, it remains an outstanding issue of whether 
the nematic order above the magnetostructural transi- 
tion is orbitally or magnetically driven. Thermodynamic 
measurement cannot distinguish them as both transitions 
fall into the Ising universality class. Indeed, increasing 
sample quality does drive the magnetic transition closer 
to the structural transitioni^ In our point of view, these 
two degrees of freedom actually coexist and cooperate 
with one another. A nonzero orbital order can certainly 
induce a spin nematicity without long-range magnetic or- 
der, and vice versa. It will be interesting if future INS 
experiments can settle the onset temperature of the spin 



nematic order and compare it to that of the orbital order 
measured by ARPESj^ However, to account for the large 
anisotropy of the magnetic exchanges orbital degrees 
of freedom have to be considered. 

One ingredient we ignored in our theory is the role 
of the structural distortions, which in principle should 
enhance the orbital order. We also neglected possible 
strong-correlation physics, which may cause an orbital- 
selective Mott transitioui^ To conclude, orbital order is 
an important part of the minimal model of iron-based 
superconductors. It induces strong exchange anisotropy 
between the local moments, and leads to an electron ne- 
matic order on the itinerant level, governing the physics 
over a very large energy scale. The study of the interplay 
between the orbital and magnetic orders should shed new 
insight into the mechanism of superconductivity in this 
multi-orbital system. 
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